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We present significant evidence of halo assembly bias for SDSS redMaPPer galaxy clusters in the 
redshift range [0.1,0.33]. By dividing the 8,648 clusters into two subsamples based on the average 
member galaxy separation from the cluster center, we first show that the two subsamples have very 
similar halo mass of M 2 oom — 1.9 x 10^'* based on the weak lensing signals at small radii 

R 10 /i“^Mpc. However, their halo bias inferred from both the large-scale weak lensing and the 
projected auto-correlation functions differs by a factor of ~1.5, which is a signature of assembly 
bias. The same bias hypothesis for the two subsamples is excluded at 2.5cr in the weak lensing and 
4.4cr in the auto-correlation data, respectively. This result could bring a significant impact on both 
galaxy evolution and precision cosmology. 


Since massive cluster-sized halos of cold dark matter 
(CDM) emerge from the rarest peaks in the primordial 
Gaussian random density field mm, their clustering am¬ 
plitudes at large scales is highly biased compared to the 
underlying mass distribution mg. In the standard ver¬ 
sion of the halo model, the clustering amplitude depends 
only on halo mass [see e.g.,[7]. 

Do secondary parameters other than the halo mass af¬ 
fect the clustering amplitude? Do observations of galaxy 
clusters reveal halo assembly bias, an effect predicted 
both by analytical theory and simulations mn]? 

In this Letter, we use a combination of the clustering 
and the weak gravitational lensing of clusters and present 
the first signihcant evidence of a difference in the large 
scale bias for cluster samples of the same mass. We divide 
our sample of galaxy clusters into two based on the aver¬ 
age projected radial separation of member galaxies, and 
investigate the difference in their clustering amplitude on 
large scales. In Ref. [13], we will show that the two clus¬ 
ter subsamples have different mass accretion rates and 
hence different assembly histories, confirming these ob¬ 
servations to be strong evidence for halo assembly bias. 
The distinguishing feature of our analysis is to use weak 
lensing to verify that the subsamples have similar halo 
masses, but different halo biases. There have been sev¬ 
eral claims for the evidence of assembly bias on galaxy 
scales [e.g., dnig. However, Ref. ng found that the dif¬ 
ference in clustering properties could all be explained as 
due to difference in halo mass or contamination by satel¬ 
lite galaxies and concluded that there was no significant 
evidence of the assembly bias for galaxy-scale halos. 

Throughout this Letter, we adopt a flat ACDM cosmo¬ 


logical model with matter density parameter Dm = 0.27 
and the Hubble parameter h = 0.7. 

Cluster subsamples - We use the publicly available cat¬ 
alog of galaxy clusters identified from the SDSS DR8 
photometric galaxy catalog by the red-sequence Matched- 
filter Probabilistic Percolation (redMaPPer) cluster find¬ 
ing algorithm [v5.10 at jSnj, also see Refs. lTHlfT^ for de¬ 
tails]. redMaPPer uses the ugriz magnitudes and their 
errors, to group spatial concentrations of red-sequence 
galaxies at similar redshifts into clusters. For each clus¬ 
ter, the catalog contains an optical richness estimate A, 
a photometric redshift estimate Zx, as well as the posi¬ 
tion and probabilities of 5 candidate central galaxies Pcen 
[IB] . A separate member galaxy catalog provides a list 
of members for each cluster, each of which is assigned a 
membership probability, Pmem [H] • 

We use a sample of clusters with 20 < A < 100 and 
0.1 < 2 ;), < 0.33. The richness cuts ensure a pure and 
statistically meaningful sample of clusters at all rich¬ 
ness bins, while the redshift cuts select a nearly volume- 
limited sample of clusters [TH|, resulting in a sample of 
8,648 clusters. For the weak lensing and clustering mea¬ 
surements, we use 100 times as many random points as 
real clusters, incorporating the survey geometry, depth 
variations, and distributions of cluster redshift and rich¬ 
ness [seeHUmi for details on the use of random catalogs]. 

As a proxy for the assembly history of the clusters, we 
use the average projected separation of member galaxies 
from the cluster center, {Rmem}- For each cluster, we 
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where the summation runs over all member galaxies, and 
Ri is the projected separation of the *-th member from 
the cluster center. Throughout this Letter use the po¬ 
sition of the most probable central galaxy in each cluster 
region as a proxy of the cluster center. We employ 14 
equally spaced bins both in redshift and A and obtain a 
spline fit for the median of (i?mem) as a function of red- 
shift and richness. We then define the two subsamples by 
the upper and lower halves of clusters in each bin of rich¬ 
ness and redshift space [23]. The ratio of (i?mem) for a 
bootstrapped realization of galaxy cluster pairs from the 
large- and the small-(i?niem) samples selected within the 
same richness and redshift has a distribution with me¬ 
dian 1.18lQ;|!)g (the errorbars correspond to the 16th and 
84th percentile). The large- and small-subsam¬ 
ples consist of 4,235 and 4,413 clusters, respectively. By 
construction, the two subsamples have almost identical 
distributions of redshift and richness. 

Weak tensing measurements - The weak gravitational 
leasing (WL) effect on the shapes of background galax¬ 
ies can be used to measure the average mass distribu¬ 
tion around galaxy clusters. We use the shape catalog of 
Ref. [ 21 ], which is based on the photometric galaxy cat¬ 
alog from the SDSS DR8 for this purpose. The galaxy 
shapes are measured by the re-Gaussianization technique 
[25], and the systematic uncertainties involved in the 
shape measurements have been investigated in great de¬ 
tail in Ref. [H] . The redshifts of source galaxies are esti¬ 
mated based on the photo-z code ZEBRA [211 jUj . The 
accuracy of the photometric redshift is not crucial for our 
study, because the populations of source galaxies used to 
compare the WL signals of our two cluster subsamples 
are identical. To measure the cluster WL profiles, we use 
the same method as described in Refs. [^1^ . 

The top and middle panels of Fig. show AE(i?), the 
excess surface mass density at a given projected radius 
R [21] for the large- and small-(i?mem) subsamples of 
clusters, respectively. The covariance matrix for each of 
the measurements was estimated based on 83 jackknife 
regions of approximately equal area covering the SDSS 
footprints (63 and 20 for the northern and southern hemi¬ 
sphere footprints, respectively). The figure shows that 
the WL signals of the two subsamples have very similar 
amplitudes at small radii, i? ^ 10 /i“^Mpc, and conse¬ 
quently similar average halo masses. However, the WL 
signals on larger scales, 10 ;< ^Mpc] < 50, display 

a significant difference, a signature of assembly bias, as 
explicitly shown in the bottom panel. 

We perform halo model fits to the measurements of the 
WL signal of each cluster subsample. Following Ref. ^5] . 
we employ a simple six parameter model fit to the WL 



FIG. 1. Halo mass consistency and assembly bias from the 
WL signal.: The data points with errorbars in the top and 
middle panels show the excess surface mass density profile as 
a function of the cluster-centric projected radius (in comoving 
units), obtained from the WL measurements for the large- and 
small-(Amem) subsamples of redMaPPer clusters (see Eq. [^, 
respectively. The points from the top panel are reproduced 
in semi-transparent color in the middle panel for compari¬ 
son. The mass profiles at small radii f? 10 /i~^Mpc appear 
to have similar amplitudes with slightly different shapes, but 
show a difference in amplitude at i? 10 h,“^Mpc, as ex¬ 
pected from assembly bias. The bold solid line shows the 
best-fit halo model, the thin solid line is the centered 1-halo 
term, the dashed line is the off-centered 1-halo term, while the 
dotted line corresponds to the stellar mass contribution from 
the central galaxy. Comparison between the dot-dashed lines 
in the two panels implies that the 2-halo term contributions, 
which arise from the average mass distribution surrounding 
the clusters, are different by a factor of 1.6. The bottom 
panel shows the ratios of the lensing signals, highlighting a 
clear deviation from unity at R 15 Mpc/h. 


signal, 

AAj(R, AfgoOm 5 t: 200 m; Qcen 5 CTofF 5 , 5) 

= <?cen^A]^^''^(i?; M200m, C200m) 

+ (1 — 9cen)AE^^'^’°®(ii; MgoOm, C200m, Ooff) 

+AE*{R; M*) -h AE2-h‘^'°(R; b). (2) 
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FIG. 2. The posterior distributions of halo model parameters 
given the WL signals for each of the two cluster subsamples 
shown in Fig. The distributions include marginalization 
over nuisance parameters which correspond to off-centering 
effects and stellar mass contribution from the central galaxy 
(see text for details). The parameters M 2 oom and C 2 oom are 
the halo mass and concentration parameters that specify the 
average NFW profile of the clusters (1-halo term), while b is 
the linear halo bias of the cluster subsample. The posterior 
distributions show the large- and small-subsamples 
have similar halo mass, but display a significant difference in 
their bias parameters. 


The first term corresponds to the halo mass profile for 
the fraction gcen of clusters whose centers have been 
correctly identified, while the second term corresponds 
to the clusters with misidentified centers. We assume 
that the halo mass profile is a smoothly truncated ver¬ 
sion of the NFW profile EDI, proposed in Ref. ED, 
specified by the halo mass and concentration parame¬ 
ter, Maoom and C 2 oom ED]- We adopt r 2 oom = 2.6 for 
the smoothing kernel, although we confirmed that our 
result of similar masses for the two subsamples is not 
sensitive to the chosen value of T 2 oom- We simply con¬ 
sider a single mass bin for host halos. We assume that 
the normalized profile of the positions of off-centered 
clusters with respect to their true center is given by 
Uos{r) oc exp[—r^/(2aQfji?|gQ^)], where Ooff describes 
the ratio of the off-centering radius to i? 200 m- We also 
truncate the off-centering profile to zero at r > i? 200 m- 
The third term, AS* = M*/(7ri?^), models a possible 
stellar mass contribution from the central galaxies assum¬ 
ing a point mass. The fourth term AE^“*'’^*°(i?) models 
the lensing contribution arising from the two-point cor- 



FIG. 3. Halo assembly bias from the projected clustering 
signal: The projected auto-correlation function of clusters for 
each of the large- and small-(Rmem) subsamples, relative to 
that of the full sample (i.e. all clusters). The clustering sig¬ 
nals at large separations, R ^ 2 /i“^Mpc, show a significant 
difference, which is consistent with the WL measurements in 
Fig.0 as shown by the shaded regions. 


relation function between the clusters and the surround¬ 
ing mass distribution. We employ the model given as 
AE^^(R) = b J{kdk/ 2 TT)p^P^lk;Zc\)J 2 ikR), where 
is the mean mass density today, b is the linear bias pa¬ 
rameter, and P^{k;Zc\) is the linear mass power spec¬ 
trum at the averaged cluster redshift Zd = 0.24, for the 
ACDM model. 

We explore the posterior distribution of the param¬ 
eters given our measurements using the Markov Chain 
Monte Carlo (MCMC) technique [53]. We use fiat pri¬ 
ors for all the parameters: M 2 oom/[ 10 ^^ h~^MQ] G 
[0.5,50], C 2 oom G [1,10], q G [0,1], Ooff G [10-4, IJ^, 
M*/[1042 G [0,10], and b G [0,10]. In Fig. § 

we show the posterior distributions of the parameters 
lW 200 m, C 2 oom as well as b, comparing results for the small- 
and large-(i?mem) subsamples, after marginalization over 
the off-centering parameters and the stellar mass contri¬ 
bution [34]. The halo masses are consistent with each 
other within the errorbars: M2oom/[10^4 /i-^M©] = 
or 1.88to;i8 for the small- and large-(i?niem) 
subsamples, respectively. The concentration parameters 
have strong degeneracies with the off-centering parame¬ 
ters, but turn out to be similar for the two subsamples 
after the marginalization. The halo bias parameters are 
b = 2.17 ± 0.31 and respectively. The ratio 

^ & 2.5cr deviation 

from unity. For comparison, even if we take the halo 
masses for the two subsamples at the extreme ends of 
their posterior distributions within their 95% C.L. inter¬ 
val, the halo bias model of Ref. [5] predicts that the ratio 
is at most 1.13. 
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Clustering measurements - We now consider the auto¬ 
correlation functions of clusters in the two subsamples to 
further confirm the difference in the large-scale bias in 
the WL signals. To avoid redshift-space distortions, we 
consider the projected correlation function 

wpiR) = 2 dne(i?,n), (3) 

Jo 

where R and If are the projected and line-of-sight sep¬ 
arations between cluster pairs, and ^(i?. If) is the three- 
dimensional correlation function. We compute the three- 
dimensional correlation function ^(i?. If) using the Tandy 
& Szalay estimator |35j , and replace the integral over line- 
of-sight by discrete summation with AH = 1 /i“^Mpc. 

In Fig. we show the ratio of the projected auto¬ 
correlation functions measured from our subsamples, rel¬ 
ative to that of the parent sample (i.e., all the clus¬ 
ters), along with jackknife error estimates. Over the 
range of separations we have considered, the two sub¬ 
samples show significantly different clustering amplitudes 
than the parent sample, giving an independent con¬ 
firmation of assembly bias. To quantify the signifi¬ 
cance, we ht a constant parameter model, Coi to C(^) = 

^^smaii-(R„,em> ^ accounting for the covari¬ 
ance between the data points. We obtain Co = 1.40±0.09 
(1.34±0.12 if we restrict to i? > 10/i“^Mpc), correspond¬ 
ing to a 4.4cr deviation from Co = 1- Using both the ratios 
of the clustering signal and the lensing signals at the three 
outermost bins in Fig. we obtain Co = 1.41 ± 0.09 with 
the significance of 4.7cr, after accounting for the cross¬ 
covariance between the clustering and lensing signals. To 
be explained by halo mass alone, the bias ratio of 1.4 
(1.6) implied from the clustering (lensing) measurements 
requires a factor of 3 (4) difference in the halo masses, 
which is disfavored by the 1-halo lensing signals. 

The shaded regions in Fig. show that the difference 
from the clustering amplitudes is consistent with the WL 
measurements, within the errorbars. The shaded regions 
were obtained with the results of fits to the lensing signals 
of the subsamples and the full sample as in Fig. They 
show the posterior distribution of the square of bias ratio 
of each subsample to the full sample, after marginaliza¬ 
tion over other parameters. The dashed lines correspond 
to the value of the maximum posterior probability. 

Conclusion and discussion — A combination of WL and 
clustering measurements of redMaPPer galaxy clusters 
show evidence of the assembly bias for cluster-sized halos 
with mass ^ 2 x 10^^ with the significance level 

of 2.5cr and 4.4 ct for each of two observables. One of the 
implications of this measurement is that the halo mass 
based on large-scale clustering amplitude is sensitive to 
the cluster selection function, which should be taken into 
account for precision cosmology. 

Could the large difference in halo bias arise from ori¬ 
entation biases in the cluster identification algorithm? 


Could the small-(i?mem) clusters preferentially reside in 
smaller mass halos and consequently have smaller bias, 
but their weak lensing mass is large due to a filament 
aligned along the line-of-sight? For the halo bias differ¬ 
ence to arise just from a halo mass difference, the sub¬ 
samples have to differ in halo mass by a factor of ^ 3 — 4, 
which is quite unlikely to occur just by projection effects 
[see e.g., [Ml EZ]- Clusters affected by projection (intrin¬ 
sically lower mass systems) are expected to cause an in¬ 
crease, contrary to a decrease in (i^mem) [3H]- We have 
explicitly verified that the line-of-sight velocity disper¬ 
sions of member galaxies around the two cluster subsam¬ 
ples are very similar as well. Furthermore, our prelimi¬ 
nary investigations indicate that a number of properties 
of central galaxies (such as their stellar masses, stacked 
spectra, velocity dispersions) in the two subsamples also 
do not show any significant differences. 

It is worth exploring whether the amplitude of the 
observed assembly bias is consistent with predictions in 
ACDM cosmologies. Ref. uni used A-body simulations 
to study the halo assembly bias in ACDM models, and 
found that the bias difference for rare objects such as our 
clusters, when subdivided by the halo concentration, is 
~ 1.25, somewhat smaller than our finding (difference of 
1.40 ± 0.09 for the clustering measurement). However, 
there are differences in our method: we used the distri¬ 
bution of member galaxies (which should correspond to 
subhalo locations) for the subdivision, the lensing mea¬ 
surements do not show a clear difference in the halo con¬ 
centration for the two subsamples (Fig.[^, and the scales 
we consider include the mildly non-linear regime. Our 
preliminary analysis using a high-resolution A-body sim¬ 
ulation, indicates qualitative agreement between the sign 
and strength of the signal in the simulations and in the 
data. These results will be presented elsewhere [T3] . 

A proper comparison between data and simulations 
will require a careful study of the cluster selection func¬ 
tion, as well as observational effects due to line-of-sight 
projections and the photo-z uncertainties, which can be 
studied by running the same cluster-finding algorithm 
on realistic simulated galaxy catalogs with full halo as¬ 
sembly histories. Additionally, on the observational side, 
one could improve the signal-to-noise ratio of the mea¬ 
surement by using the cross-correlations of the cluster 
samples of this study with a sample of galaxies such as a 
catalog of luminous red galaxies that has a much higher 
spatial number density. We defer these topics to a study 
in the near future. 
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